SIMU=0;CD_REAL_20_0=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_0[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_0=c(CD_REAL_20_0,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_0,'CD_REAL_20_0.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_1=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_1[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_1=c(CD_REAL_20_1,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_1,'CD_REAL_20_1.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_2=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_2[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_2=c(CD_REAL_20_2,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_2,'CD_REAL_20_2.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_3=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_3[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_3=c(CD_REAL_20_3,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_3,'CD_REAL_20_3.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_4=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_4[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_4=c(CD_REAL_20_4,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_4,'CD_REAL_20_4.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_5=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_5[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_5=c(CD_REAL_20_5,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_5,'CD_REAL_20_5.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_6=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_6[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_6=c(CD_REAL_20_6,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_6,'CD_REAL_20_6.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_7=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_7[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_7=c(CD_REAL_20_7,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_7,'CD_REAL_20_7.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_8=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_8[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_8=c(CD_REAL_20_8,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_8,'CD_REAL_20_8.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_9=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_9[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_9=c(CD_REAL_20_9,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_9,'CD_REAL_20_9.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_10=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_10[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_10=c(CD_REAL_20_10,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_10,'CD_REAL_20_10.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_15=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_15[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_15=c(CD_REAL_20_15,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_15,'CD_REAL_20_15.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_20=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_20[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_20=c(CD_REAL_20_20,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_20,'CD_REAL_20_20.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_25=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_25[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_25=c(CD_REAL_20_25,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_25,'CD_REAL_20_25.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_30=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_30[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_30=c(CD_REAL_20_30,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_30,'CD_REAL_20_30.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_35=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_35[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_35=c(CD_REAL_20_35,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_35,'CD_REAL_20_35.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_40=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_40[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_40=c(CD_REAL_20_40,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_40,'CD_REAL_20_40.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_45=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_45[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_45=c(CD_REAL_20_45,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_45,'CD_REAL_20_45.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_50=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_50[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_50=c(CD_REAL_20_50,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_50,'CD_REAL_20_50.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_55=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_55[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_55=c(CD_REAL_20_55,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_55,'CD_REAL_20_55.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_60=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_60[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_60=c(CD_REAL_20_60,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_60,'CD_REAL_20_60.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_65=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_65[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_65=c(CD_REAL_20_65,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_65,'CD_REAL_20_65.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_70=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_70[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_70=c(CD_REAL_20_70,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_70,'CD_REAL_20_70.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_75=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_75[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_75=c(CD_REAL_20_75,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_75,'CD_REAL_20_75.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_80=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_80[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_80=c(CD_REAL_20_80,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_80,'CD_REAL_20_80.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_85=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_85[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_85=c(CD_REAL_20_85,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_85,'CD_REAL_20_85.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_90=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_90[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_90=c(CD_REAL_20_90,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_90,'CD_REAL_20_90.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_95=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_95[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_95=c(CD_REAL_20_95,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_95,'CD_REAL_20_95.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_100=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_100[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_100=c(CD_REAL_20_100,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_100,'CD_REAL_20_100.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_110=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_110[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_110=c(CD_REAL_20_110,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_110,'CD_REAL_20_110.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_120=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_120[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_120=c(CD_REAL_20_120,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_120,'CD_REAL_20_120.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_130=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_130[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_130=c(CD_REAL_20_130,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_130,'CD_REAL_20_130.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_140=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_140[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_140=c(CD_REAL_20_140,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_140,'CD_REAL_20_140.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_150=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_150[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_150=c(CD_REAL_20_150,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_150,'CD_REAL_20_150.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_160=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_160[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_160=c(CD_REAL_20_160,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_160,'CD_REAL_20_160.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_170=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_170[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_170=c(CD_REAL_20_170,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_170,'CD_REAL_20_170.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_180=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_180[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_180=c(CD_REAL_20_180,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_180,'CD_REAL_20_180.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_190=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_190[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_190=c(CD_REAL_20_190,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_190,'CD_REAL_20_190.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_200=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_200[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_200=c(CD_REAL_20_200,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_200,'CD_REAL_20_200.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_210=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_210[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_210=c(CD_REAL_20_210,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_210,'CD_REAL_20_210.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_220=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_220[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_220=c(CD_REAL_20_220,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_220,'CD_REAL_20_220.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_230=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_230[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_230=c(CD_REAL_20_230,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_230,'CD_REAL_20_230.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_240=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_240[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_240=c(CD_REAL_20_240,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_240,'CD_REAL_20_240.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_250=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_250[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_250=c(CD_REAL_20_250,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_250,'CD_REAL_20_250.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_260=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_260[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_260=c(CD_REAL_20_260,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_260,'CD_REAL_20_260.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_270=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_270[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_270=c(CD_REAL_20_270,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_270,'CD_REAL_20_270.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_280=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_280[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_280=c(CD_REAL_20_280,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_280,'CD_REAL_20_280.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_290=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_290[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_290=c(CD_REAL_20_290,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_290,'CD_REAL_20_290.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_300=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_300[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_300=c(CD_REAL_20_300,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_300,'CD_REAL_20_300.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_310=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_310[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_310=c(CD_REAL_20_310,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_310,'CD_REAL_20_310.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_320=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_320[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_320=c(CD_REAL_20_320,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_320,'CD_REAL_20_320.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_330=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_330[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_330=c(CD_REAL_20_330,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_330,'CD_REAL_20_330.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_340=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_340[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_340=c(CD_REAL_20_340,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_340,'CD_REAL_20_340.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_350=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_350[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_350=c(CD_REAL_20_350,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_350,'CD_REAL_20_350.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_360=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_360[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_360=c(CD_REAL_20_360,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_360,'CD_REAL_20_360.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_370=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_370[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_370=c(CD_REAL_20_370,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_370,'CD_REAL_20_370.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_380=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_380[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_380=c(CD_REAL_20_380,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_380,'CD_REAL_20_380.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_390=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_390[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_390=c(CD_REAL_20_390,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_390,'CD_REAL_20_390.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_400=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_400[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_400=c(CD_REAL_20_400,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_400,'CD_REAL_20_400.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_410=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_410[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_410=c(CD_REAL_20_410,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_410,'CD_REAL_20_410.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_420=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_420[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_420=c(CD_REAL_20_420,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_420,'CD_REAL_20_420.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_430=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_430[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_430=c(CD_REAL_20_430,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_430,'CD_REAL_20_430.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_440=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_440[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_440=c(CD_REAL_20_440,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_440,'CD_REAL_20_440.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_450=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_450[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_450=c(CD_REAL_20_450,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_450,'CD_REAL_20_450.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_460=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_460[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_460=c(CD_REAL_20_460,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_460,'CD_REAL_20_460.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_470=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_470[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_470=c(CD_REAL_20_470,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_470,'CD_REAL_20_470.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_480=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_480[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_480=c(CD_REAL_20_480,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_480,'CD_REAL_20_480.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_490=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_490[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_490=c(CD_REAL_20_490,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_490,'CD_REAL_20_490.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_20_500=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_500[,w];y2=PROFILES_G2_20[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_20_500=c(CD_REAL_20_500,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_20_500,'CD_REAL_20_500.txt',row.names=FALSE,col.names=FALSE)

